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Abstract. We study the problem of constructing systems of hyperbolic conservation laws in one 
space dimension with prescribed eigencurves, i.e. the eigenvector fields of the Jacobian of the flux are 
given. We formulate this as a typically overdetermined system of equations for the eigenvalues-to-be. 
Equivalent formulations in terms of differential and algebraic-differential equations are considered. 
The resulting equations are then analyzed using appropriate integrability theorems (Frobenius, 
Darboux and Cartan-Kahler). We give a complete analysis of the possible scenarios, including 
examples, for systems of three equations. As an application we characterize conservative systems 
with the same eigencurves as the Euler system for f-dimensional compressible gas dynamics. The 
case of general rich systems of any size (i.e. when the given eigenvector fields are pairwise in 
involution; this includes all systems of two equations) is completely resolved and we consider various 
examples in this class. 
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1. Introduction 

Consider a system of n conservation laws in one space dimension written in canonical form, 
(1.1) u t + f(u) x = 0. 

Here the unknown state u = u(t, x) £ W 1 is assumed to vary over some open subset C W 1 and 
the flux / = (f 1 , . . . , f n ) T is a nonlinear map from f2 into K n . The eigenvalues and eigenvectors of 
the Jacobian matrix Df(u) provide information that is used to solve the Cauchy problem for (jl.lj) . 
In particular, the geometric properties of the integral curves of eigenvector fields of Df play a key 
role. Together with the so-called Hugoniot curves (see below for definition) these form wave curves 
in li-space that are used to build solutions to (|l.ip . 

There is a well developed theory for the Cauchy problem for a large class of systems (jl.ip in 
the near equilibrium regime u(0, x) ~ constant. Glimm [T7J established global-in-time existence of 
weak solutions for data with sufficiently small total variation. By now there are several methods 
available to obtain these solutions: Glimm's original random choice method [17], Liu's deterministic 
version [29] , the wave- front tracking schemes of DiPerna [2] , Bressan [5] and Risebro [33] , and the 
wave tracing methods of Bianchini and Bressan [1] , [2] , [3] . 

Single equations and systems of two equations enjoy particular properties that yield global exis- 
tence beyond the perturbative regime, [26], [IB], [30], pD], [30], pj)], [8]. On the other hand there is 
no general existence result available for "large" data when the system has three or more equations 
(e.g. the Euler system for compressible gas dynamics). Indeed, examples of blowup in total varia- 
tion and/or sup norm are known even for genuinely nonlinear and strictly hyperbolic systems, 
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[21) . |22j . pQ. These examples show that the global behavior of the wave curves in u-space is 
related to finite time blow up, [H], [21], [22], [I]. It is of interest to have more examples of this 
type. In particular one would like to know if a physical system can display similar behavioiQ- In 
searching for systems whose wave curves have special properties one is naturally lead to ask what 
freedom one has in prescribing such curves. 

In this paper we consider the situation where we are given a frame of n linearly independent 
vector fields and their integral curves. We then ask if there are any systems of the form (|l.ip with 
the property that the given vector fields are the eigenvector fields of Df(u). 

It turns out that there is always a one-parameter family of trivial solutions, and that there are 
cases where no non-trivial solutions exist. More generally we are interested in knowing "how many" 
solutions there are, in terms of how many constants and functions appear in a general solution. The 
basic tool for this sort of questions is provided by the integrability theorems for overdetermined 
systems of PDEs such as Frobenius, Darboux and Cartan-Kahler theorems. The latter requires 
re- writing the system of PDE's as an exterior differential systems (EDSs) [7], [20] . 

The case of 2 x 2-systems was treated in PjQ; it is also covered by the analysis of rich systems 
in Section HJ On the other hand we shall see that already for systems with three equations there 
are several possible scenarios. Before giving a precise formulation of the problem we review some 
relevant background material. 

Notation: We denote the (i,j')-entry (i.e., the element in the ith row and the jth column) of 
an m x n-matrix A by A 1 -. Superscript T denotes transpose. Summation convention is not used. 

1.1. Conservation laws in one space dimension. We consider hyperbolic systems of conserva- 
tion laws in one spatial dimension (11. ip . i.e. the Jacobian matrix Df{u) is diagonalizable over M. at 
each state u G VL. The system is strictly hyperbolic in f2 provided the eigenvalues X l (u) of Df(u) 
are real and distinct: 

(1.2) AV) < ••• < X n (u), Wuefl. 

Let us for now fix a choice of the associated right and left eigenvectors Ri(u) and L l (u) of Df(u). 
These are considered as column and row vectors of functions, respectively, and we write 

Ri(u)= [Rj («),..., R? (u)] T , V (u) = [L\ («),..., 4(n)] . 

We refer to the Ri{u) as the eigenfields and their integral curves in u-space as eigencurves. Diago- 
nalizing Df we have 

(1.3) Df(u) = R(u) A(u) L(u) , 
where 

R(u) = [i?i(u) | ••• | Rn(u)} , A(«) = diagfA 1 ^) . . . X n (u)] , 

and 

" LHu) " 



L(u) = R(uY 



L n (u) 

Note that in setting L(u) = i?(u) _1 we have introduced the normalization Ri(u) • ]J(u) = b\ 
(Kronecker delta). 

Next consider the initial value problem for (jl.ip where the data consist of two constant states 
separated by a jump at x = 0, 

u- , x < 
u+ , x > . 



(1.4) u (x) 



^This refers to the strictly hyperbolic regime. It is well-known that failure of strict hyperbolicity can cause singular 
behavior (blowup in total variation or L°°), even in physical systems. See [41], |37| and references therein. 
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This is the so-called Riemann problem and its solution serves as a building block for more general 
solutions, [27], [T7J. For a sufficiently smooth flux / whose characteristic fields are genuinely 
nonlinear or linearly degenerate in the sense of Lax [27], it is well known that through every 
strictly hyperbolic state u there exist n locally defined and C 2 smooth wave curves. These curves 
collectively provide self-similar solutions to Riemann problems; for details see [6], [12], [38]. Each 
wave curve is locally made up of two components with second order contact at the base point u: 
the rarefaction states that are part of the eigencurves, and the shock states that are part of the 
Hugoniot locus {u £ Q \ 3 s E R : f (u) — f (u) = s(u — u) } . The geometry of these curves in state 
space thus provide information about the solutions to (|1.1|) . 

1.2. Connections on frame bundles. Given an n-dimensional smooth manifold M we let X(M) 
and X*{M) denote the set of smooth vector fields and differential 1-forms on M, respectively. A 
frame {n, . . . , r n } is a set of vector fields which span the tangent space T p M at each point p € M. 
A coframe {I 1 , . . . ,£ n } is a set of n differential 1-forms which span the cotangent space T*M at 
each point p £ M. The coframe and frame are dual if £ l {rj) = <5j (Kronecker delta). If u 1 , . . . ,u n 
are local coordinate functions on M, then {^r, • • • , gfrr} is the corresponding local coordinate 
frame, while {du 1 , . . . , du n } is the dual local coordinate coframe. For a given frame {ri, . . . , r n } the 
structure coefficients c^. are defined through 

n 

(1-5) [n,r j ]=J2 c ij r k, 

k=l 

and the dual coframe has related structure equations given by 
(1.6) d£ k = -^c%t A£ j . 

i<j 

It can be shown that there exist coordinate functions w l , . . . ,w n on f2 such that = ^j-, i = 
1 . . . ,n, if and only if r±,...,r n commute, i.e. all structure coefficients are zero. Next, an affine 
connection V on M is an M-bilinear map 

X(M) x X(M) -> X(M) (X, Y) i-> V X Y 

such that for any smooth function / on M 

(1-7) V fX Y = fV x Y, V x (fY) = (Xf)Y + fV x Y. 

By R-bilinearity and (jl.7p the connection is uniquely defined by prescribing it on a frame: 

n 
k=l 

where the smooth coefficients T^- are called connection components, or Christoffel symbols, relative 
to the frame {r l5 . . . ,r n }. Any choice of a frame and n 3 functions T^-, i,j,k = 1, . . . ,n, defines 
an affine connection on M. A change of frame induces a change of the connection components, 
and this change is not tensorial. E.g., a connection with zero components relative to a coordinate 
frame, may have non-zero components relative to a non-coordinate frame. 

Given a frame {n, . . . , r n } with associated Christoffel symbols T^- and dual frame {I 1 , . . . , £ n }, 

we define the connection 1-forms fM- by 

(i-s) /4==ErJA 

k=l 
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In turn, these are used to define two important tensor- fields: the torsion 2-forms 

n 

(1.9) T*:=df +^ / 4A£ fe = J^Ar, i = l,...,n, 

A;=l k<m 

and the curvature 2-forms 

(1.10) Rf := ^" + £ 4 A ^ = £ i4„/ A f» . 

fc=l k<m 

Here 

/ 1 i i \ mi -pi _ -p i _ i 

V / fcm km mk ^km 

n 

(i-i2) flL, = a(ry-r m (r{J+^(ri s r^-r^-4 m r^ 

are components of the torsion and curvature tensors respectively, and these do change tensorially 
under a change of frame. We can write equations (jl.9p and (jl.lOp in the compact matrix form 

(1.13) T = d£ + fiA£, R = dfi + \i A n 

where £ = (£ 1 , . . . ,£ n ) T , T = (T 1 , . . . , T n ) T , and R and fi are the matrices with components R^ 
and \i\ respectively. The connection is called symmetric if the torsion form is identically zero and 
it is called flat if the curvature form is identically zero. Equivalently: 

(1.14) dt = —pL A £ (Symmetry), d/j, = —pi A [i (Flatness). 
In terms of Christoffel symbols and structure coefficients this is equivalent to 

(1.15) 4 m = T km- r mk (Symmetry) 
and 

n 

(1.16) r m (T j ki ) - r k {P mi ) = ]T (r^ - Ti s T s ki - c s km P m ) (Flatness). 

s=l 

One can also show that a connection V is symmetric and flat if and only if in a neighborhood of 
each point there exist coordinate functions u \ . . . , u n with the property that the Christoffel symbols 
relative to the coordinate frame are zero: V_§_t£j = for all i,j = 1, . . . ,n. 

1.3. Formulation of the problem. Returning to (jl.ip we next provide a precise statement of 
the "inverse" problem of constructing flux functions / whose geometric properties are given. There 
are various ways to formulate such problems. One could prescribe families of curves which are then 
required to be the Hugoniot loci, the eigencurves, or the wave curves for a system of conservation 
laws (jl.ip . One might also consider giving combinations of these. The most direct formulation is 
obtained by prescribing the eigenfields (equivalently, the eigencurves), and this is what we do here. 

We will be working locally near a fixed base point u in an open set O C M n , with Q smoothly 
contractible to a point. Throughout u = (u , . . . ,u n ) will denote a fixed system of coordinates on 
a neighborhood of u. We assume that we are given n linearly independent column n- vectors Ri(u) 
(the eigenfields to be) on ft, and we define 
(1.17) 

R{u) := [Riiu) | ••• | R n (u)} = {R){u)) itJ , L{u) := tf(u)- 1 = 



L\u) 



L n (u) 



(£}(«)) 



The problem we consider may be formulated as follows: 
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Problem 1. Consider an open subset O C W 1 equipped with a coordinate system u = (u , . . . ,u n ). 
Given a point u £ together with a (sufficiently smooth) frame {Ri(u), . . . , R n (u)} on fi. 

Then: find n real functions A 1 (u), . . . , X n (u) defined on a neighborhood Li C U of u such that, 
with A(u) := diag[A 1 (u) . . . X n (u)], the matrix 



is the Jacobian matrix with respect to u of some map f :U — > ffi n . We are further interested in how 
large the set of solutions is, i.e. how many arbitrary constants and functions appear in a general 
solution A (it), ... , X n (u). 

A solution X(u) = (A 1 (it), . . . , X n (u)) to Problem Q] generates a flux / and an associated system 
of conservation laws (jl.ip in which u , . . . , u n are the conserved quantities. We note that we do not 
impose strict hyperbolicity for solutions to our problem; indeed, we shall see that there are cases 
where two or more eigenvalues X 1 must necessarily coalesce. 

Below we derive a system of partial differential equations (PDEs) for the eigenvalues A*, where the 
coefficients are expressed in terms of the components of Ri and L 1 . The system may be formulated 
in different ways, all of which we refer to as the X-system. Typically this will be an overdetermined 
system of linear, variable coefficients PDEs. It will turn out that the most useful formulation of 
the A-system is as an algebraic- differential system; see Section [2l 

Remark 1.1. Let us clarify the coordinate dependence in Problem^ The property of being a 
Jacobian with respect to a system of coordinates is not invariant under most changes of coordi- 
nates. For this reason we need to fix the coordinates u % at the outset. Of course, having started 
with a particular choice of coordinates, it may be that the resulting X-system takes a simpler form 
when formulated in other coordinates. This is the case for so-called rich systems where the "good" 
coordinates are the Riemann invariants (see Definition^ below and Section^. 

It may also be that a change of coordinates happens to preserve Jacobians. This occurs for the 
Euler system of compressible flow when written in Lagrangian variables (Example \5.1\) . In this case 
the nonlinear change of dependent variables (specific volume, velocity, total energy) (specific 
volume, velocity, entropy) preserves the conservative form of the system. As far as Problem^ is 
concerned the two forms are equivalent. 

Remark 1.2. Besides its intrinsic interest Problem{l\is motivated by possible applications where 
one searches for systems M.l\) with particular geometric properties, and associated properties (e.g. 
blowup) of solutions to Cauchy problems. In our setup we thus consider the "geometry" (the eigen- 
fields) as given, and we look for the scalar eigenvalue fields. Alternatively one could prescribe the 
eigenvalues and search for eigenfields Ri(u) that renders A(u) in U.18\) a Jacobian matrix. 

A special class of systems is given by so-called rich systems. These are systems equipped with 
a coordinate system of Riemann invariants. For definitions, and the fact that richness can be 
expressed in terms of the eigenfields, we refer to [35], and Section 7.3 in [12]. Starting from a given 
frame it is convenient to make the following, equivalent, definition: 

Definition 1. The frame {R\ (it), . . . , R n (u)} of linearly independent vector fields is said to be rich 
if each pair of vector fields is in involution: 



While it is natural to consider the class of rich systems in connection with Problem [T] we shall 
see that richness does not imply any uniformity with respect to how many constants and functions 
are needed to specify a general solution. 

Remark 1.3. The class of symmetrizable systems plays a central role in the general theory of 
conservation laws, [16], [19]. These are the systems that admit a convex "entropy" r] : — > R ; 
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(1.19) 



[Ri,Rj] G span{i?j, Rj} 



for all 1 < i, j < n . 



with an associated entropy-flux q : £1 — > R, such that an additional conservation law 

7](u) t + q(u) x = 

holds whenever u is a smooth solution of Equivalently (at least in the strictly hyperbolic 

regime), the Hessian D 2 n should define an inner product with respect to which the given eigenvec- 
tors R{ are orthogonal (Rj D 2 w Ri = for i ^ j; see [12] for details). Differently from richness, 
symmetrizability is not expressed directly in terms of the eigenfields alone. Problem [7] for sym- 
metrizable systems will be taken up elsewhere. 

1.4. Related works. Problem Q] was addressed by Dafermos [TT] for 2 x 2-systems in several space 
dimensions under the requirement that the Jacobians in the various spatial directions commute. 
Commutativity implies that the Jacobians have the same eigenfields. In [11] it was shown how 
to construct such systems for any pair of linearly independent vector fields. The case of one- 
dimensional 2 x 2-systems is covered by the analysis of rich systems in Section [Tj 

In his geometric analysis of systems of conservation laws Sevennec [36] provides a characterization 
of those quasilinear systems 

v t + A(v)v x = 

that can be transformed to conservative form (|l.lj) by a suitable change of variables u = <f)(v). In 
particular, the characterization involves a version of what we refer to as the A-system. 

The class of rich systems has been studied by many authors. Conlon and Liu [JJ)] considered 
rich systems in connection with entropy criteria and showed that such systems are endowed with 
large families of entropies. From a different perspective the same class of systems were studied 
by Tsarev [421143]. Sevennec [36] showed that the eigenvalues of strictly hyperbolic, rich systems 
must satisfy certain restrictive conditions (see Proposition 5 in [36]). Serre [35] has performed 
a comprehensive analysis of rich systems, including building of entropies, commuting families of 
systems, and construction of rich systems. In Section 14.21 we analyze rich systems in relation to 
Problem [TJ 

1.5. Outline and summary of results. It turns out that the complete solution of Problem Q] for 
arbitrary n is quite complicated. In this paper we provide a complete solution for the case n = 3, 
as well as for rich systems for any n. (This latter class covers the case n = 2.) We also give a list 
of examples that illustrate the various types of solutions. 

In Section [2] we begin by noting some properties the A-system that follow directly from the 
formulation of Problem [TJ We then formulate three equivalent versions of the A-system, including 
an algebraic-differential system, and we record the extreme cases with minimal and maximal number 
of algebraic constraints. A few general facts are collected in Proposition 12.21 

Our solution of the n = 3 case in Section [3] reveals that the solution set of the A-system depends 
on the number of independent algebraic equations, as well as on the number of A* that appear in 
these equations. For n > 4 there seems to be no easy way to analyze completely the resulting 
cases. We provide a complete breakdown of the possible scenarios when n = 3. The algebraic part 
of the A-system now contains two, one, or no independent algebraic equations. The two extreme 
cases fall into either the trivial or rich categories. In the case of one algebraic equation the size 
of the solution set depends further on the number of A* involved in the algebraic constraint. The 
Frobenius integrability theorem can be applied when all three A* appear in the algebraic equation, 
and in this case the general solution depends on two constants. The only other possibility is that 
exactly two A* occur in the algebraic equation. This case may be analyzed by using the Cartan- 
Kahler integrability theorem, and the general solution now depends on one function of one variable 
and one constant. 

We solve Problem [TJ for the class of rich systems of any dimension in Section [H The subclass of 
rich systems without algebraic constraints are analyzed in Section f4. 21 while Section f4.3l treats the 
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more involved case of rich systems where the eigenvalues are related algebraically. In both cases 
the A-system can be analyzed using an integrability theorem of Darboux (Theorem 14. ip . 

For completeness we include statements of the various integrability theorems we apply. Concern- 
ing smoothness assumptions of the given frame we need to require analyticity when we apply the 
Cartan-Kahler theorem. On the other hand, for the cases that use the theorems of Frobenius and 
Darboux we only need to require C 2 smoothness of the given frame. 

Section [5] collects several examples that are of interest in themselves or that illustrate the different 
cases treated in Section [3] and Section UJ We start by considering the various solutions of Problem 
[Tj for the case where the given eigenfields are those of the Euler system describing one-dimensional 
compressible fluid flow (Example 15. ip . Depending on the prescribed pressure function these eigen- 
fields may form either a rich or a non-rich frame. We proceed with several more examples of 
non-rich frames on R 3 illustrating various scenarios treated in Section [3] and an example of a frame 
on K 4 with only trivial solutions for the A-system. We then give a set of examples of rich frames 
whose A-systems do not impose algebraic constraints on eigenfunctions. This includes a frame on 
M 2 , a rich orthogonal frame on M 3 and a constant frame on W 1 . We finally give two examples of 
rich frames on M 3 whose A-systems impose certain algebraic constraints on eigenfunctions. We have 
written a set of procedure^ in Maple to obtain and analyze the A-system for a given frame, and 
used it to construct some of the above examples. 

Finally, our analysis indicates that a solution of Problem [T] for general systems with n > 4 is 
rather involved, with a large number of different subcases. 

2. The A-system 

2.1. Trivial solutions and scalings. Problem [Tj always has a one-parameter family of solutions 
given by 

A x (n) = ••• = \ n (u) = A, 
where A is any real constant. We refer to these as trivial solutions. The resulting system (jl.ip 
is linearly degenerate in all families (i.e. VAj(u) • Ri{u) = 0) and any map f(u) = An + u, where 
u E M n , is a corresponding flux. While not of interest themselves, the trivial solutions show that 
any compatibility condition associated with the A-system will not rule out existence of solutions 
altogether. A trivial solution can be added to any solution of the A-system to give another solution 
of the same system. For later reference we record the following related result: 

Proposition 2.1. If A 1 (it) = • • • = X n (u) is a solution to Problem^ then their common value is 
a constant. 

Proof. If X l (u) = ■ ■ ■ = \ n (u) = \{u) then A(u) = A(u) = X(u)I nxn . For this to be a Jacobian of 
a map / : U — > W 1 we must have that 

dfi 

Vi,j = l,...,n : -^j(u) = 5ij\(u) , 

whence X(u) is a function of u l alone for each i, and thus constant. □ 

In formulating Problem [Tj we may use any (non- vanishing) scalings of the eigenfields Ri(u). That 
is, given smooth functions a 3 : — > M. \ {0}, j = l,...,n, we may set Rj(u) = a 3 (u)Rj(u), 
together with the inversely scaled left eigenfields Lj(u) := a 3 (u) 1 Lj(u). Letting R := Ra, L := 
a^ 1 L, where a{u) = diag[ai(u) . . . a n (u)], we get that R(u)A(u)L(u) is a Jacobian if and only if 
R(u)A(u)L(u) is a Jacobian. E.g., in the case of rich systems the simplest form of the A-system is 
obtained by using the versions that makes the matrix L a Jacobian with respect to the u variable 
(see Section H]). 



MAPLE code is available at http://www.math.ncsu.edu/~iakogan/mapleHTML/lambda-system.html 
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2.2. Formulating the A-system. 

2.2.1. Direct formulation. We have that a matrix A(u) = (Aj(u))ij, is a Jacobian with respect to 
it-coordinates if and only if 

(2.1) d k A)(u) = djA\(u) for all i,j,k = l,...,n with j < k, 

where d\ denotes partial differentiation with respect to U{. In the case when A(u) is given by (j!.18j) 
we set 

C l m j{u) := R l m (u)LY(u) (no summation), 

and (|2.ip may be written 

n 

(2.2) E \c l mj d k \ m - C l mk dj\ m + X m (d k C l mj - djC l mk ) = , 

m=l 

where k E {l,...,n} with j < k. This version of the A-system provides a homogeneous, 

variable coefficient system of " h near PDEs for n unknowns. For n > 3 it is thus typically 

an overdetermined system of PDEs. 

2.2.2. Formulations using 1 -forms. A simpler formulation of the A-system is obtained by using 
differential forms to express the condition that A(u) in Problem Q] is a Jacobian. By Poincare's 
lemma (recall that O is assumed smoothly contractible to a point) we have 

(2.3) A(u) is a Jacobian matrix w.r.t. u dA{u) Adu = , 

where the d-operator is applied component- wise. Applying the product rule, condition (|2.3p is thus 
equivalent to 

(2.4) {L(dR)A + dA — AL(dR)} A Ldu = , 
or 

(2.5) (dA) A {Ldu) = {A(LdR) - (LdR)A} A (Ldu) , 

where we have used that L = R . (Clearly (|2.5p is satisfied if dA = A(LdR) — (LdR)A; however, 

the associated solutions are exactly the trivial solutions (Ai, . . . , X n ) = (A, . . . , A), A G K.) 

The system (|2.5p is an equation for n-vectors of 2-forms. We proceed to write out the system in 

_d Q_- 

du 1 ' dui ■ 



u-coordinates by applying (12. 5|) to the pair of vector fields (A, ~yj)- A direct calculation yields: 



(2.6) 4(^ Ai ) " L l( d 3 Xi ) = Yl B lj\ Xm ~ Xi ) > for each * and for all 1 < / < j < n, 

where B™ := {L^ffiL 1 ) — Lf l (djL 1 )^ • R m . Not surprisingly this is again a homogeneous system 

of n ^ ^ linear PDEs for n unknowns. Thus, in this formulation of the A-system the equations 
involving derivatives of A* are such that, first, no other eigenvalue appears differentiated, and 
second, at most two derivatives of A* occur in each equation. 

To obtain a formulation of the A-system which explicitly records algebraic constraints on the 
eigenvalues A 1 , we introduce the following notation. Whereas we treat Ri(u) and L l (u) as arrays 
of functions we let ri{u), £ l (u) denote the corresponding vector fields (differential operators) and 
differential 1-forms: 

n Q n 

(2.7) ri(u) := E^(«)^[, *(«) := E 4»^ m | u ■ 

m=l m=l 



Since R(u) is assumed invertible on the vector fields ri(u), i = 1, . . . ,n, provide a frame on Q, 
and the 1-forms £ l (u),i = 1, . . . , re, provide the dual coframe on O (see Section fl . 2 1) . We define the 
n-vector of 1-forms £ by 



£ := '■ = Ldu . 

and introduce the following coefficients 

(2.8) Tfj := L k (DRj)Ri , 

where D denotes Jacobian with respect to u. A direct computation shows that they are, in fact, the 
Christoffel symbols (connection components) of the flat and symmetric connection V_a_^y = 

computed relative to the frame {n, • • • , r n }. In other words, the covariant derivatives of the frame 
vector fields are V ri rj = Y!k=i^ij r k- We note that the symmetry condition (|1.15p implies 

n n 

(2-9) h, rj ] = 4 n = £(r* - r%) n . 

k=l k=l 

This shows that if {r±, . . . , r n } is rich, then T^- = whenever k ^ {i, j}. 

Another direct calculation shows that the (fe,j)-entry of the matrix fj, := R dR = LdR of 
1-forms is given by the connection forms 



(2.10) M *:=(Ldfl)J = £r. 



k pi 
if 

i=l 

Thus (|2.5p reads 

(2.11) (dA) A£ = {Api- fiA} At. 

Again (|2.1ip is an equation of n-vectors of 2-forms. Applying each component to pairs of frame 
vector fields (rj,r\,-), i,j = 1, . . . , n, we obtain an equivalent form of the A-system as a differential- 
algebraic system: 

(2.12) ri(A') = rj,iA' A ' ) for i / j, 

(2.13) (A 4 -A fc )r^ = (A'-A*)r* fari<j, i ^ k, j ^ k, 

where there are no summations. (I2.12|) gives n(n — 1) linear, homogeneous PDEs, while (|2.13l) 
gives "( n 2 ) algebraic relations. We observe that (12. 13H is symmetric in i and j, and that all 
coefficients Ffj with i ^ j ^ k ^ i appear in (|2.13p . This form of the A-system explicitly records 
the algebraic relations inherent in (|2.2p and (|2.6p . 

Remark 2.1. The equations \2.12]) - l2.13]) appear in Sevennec's characterization of quasilinear 
systems that admit a conservative form, see [36] . For a given quasilinear system 



v t + A(v)v x = 0, A(v)eR nxn , 

Sevennec shows that there is a coordinate system in which the system is conservative if and only 
if there exists a flat and symmetric affine connection V such that its Christoffel symbols and the 
eigenvalues of A(u) satisfy 12.12]) - [2.13]) . The same system appears in Tsarev 

Remark 2.2. Note that condition \2.11]) is equivalent to 

(2.14) 6 := (dA-A/x + /xA) M = 0, 
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where is a differential form on a 2n- dimensional manifold with coordinates u , . . . , u n , A , . . . , A™. 
Moreover, the ideal I algebraically generated by in the ring of differential forms is a differential 
ideal, i. e. duj £ Z for all to £ 1. Indeed, from the flatness and symmetry properties of the connection: 

dfj, = —\x A // and dl = —\i A I , 

it follows that 

dO = ( - dA A // - A(dju) + (d//)A - fi A dA) A £ - (dA - Afi + fiA) A dl 

= ( - dA A ijl + AO A fi) - (p A /i)A - fi A dA) A £ + (dA - A^ + /xA) A/iAf 

= ( - (At A /i)A — fi A dA) A£ + (fiA) A fi A £ 

= nA(-/j>A-dA + A[i)A£ = -fiA@ = mod . 

Solving the X-systems is therefore equivalent to finding n-dimensional integral sub-manifolds X 1 = 
F l (u\...,u n ),i = l,...,n, of I. 

2.3. Observations about the rank of the algebraic sub-system ( 12.131) . We next consider 
the algebraic constraints (|2.13p . which is a system of n ( n ~ 1 K ?1 ~ 2 ) li near equations. Choosing the 
variables to be the differences 

x fc :=A fc -A 1 , k = 2,...,n, 

we rewrite (I2.13P in matrix form as 
(2.15) Nx = 0. 

Here x is the (n — l)-vector (x 2 , . . . ,x n ) T and N is a certain n ( n 1 )i n 2 ) x ( n _ i)- ma trix whose 
entries are given in terms of the T^. It is easily checked that each entry of N is either zero, a single 
Christoffel symbol (±r^), or a difference of such (r^- — T^). Furthermore, N = if and only if 
T^j = for all choices i j ^ k i. 

The number of independent algebraic constraints is given by rank(A r ). It is convenient to use 
this as a first, rough classification. However, rank(A^) does not characterize the solutions to the 
A-system (in terms of number of constants and functions in the general solution). This is clear 
already from the cases of maximal and minimal rank. We briefly consider these extreme cases: 

• rank(AQ = 0. In this case N = and there are no algebraic constraints imposed on the 
eigenvalues. This occurs if and only if T^- = for all choices of i 7^ j 7^ k 7^ i. By (12. 9j) this 
implies that [rj,rj] € spanjrj, r,}, i.e. we are in the rich case: 

rank(A r ) =0 {ri, . . . , r n } is rich . 

The case rank(A r ) = will be analyzed in more detail in Section 14.21 We will show that 
any frame on R 2 , any constant frame on M n , as well as any rich orthogonal frame on W 1 , 
falls in this category. Corresponding examples are given in Section 15.31 

The question of whether richness implies rank(A^) = is somewhat subtle. Namely, if 
we were given a rich and strictly hyperbolic system (jl.ip . then indeed rank(A r ) = 0. (See 
[T2] p. 185 for the proof that Y\- = L k (DRj)Ri = for all choices of i ^ j 7^ k 7^ i in 
this case.) However, with Problem [T] we are starting from given vector fields rj without 
insisting on strict hypebolicity. Example 15.91 and Example 15.101 show that it is possible to 
prescribe a collection of vector fields that form a rich family, without rank(A^) being zero. 
Furthermore, Example 15.91 shows that the associated A-system for frames of this type may 
admit nontrivial solutions. Thus: 

{n, . . . , r n } is rich rank(A^) = . 
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In Section 14.31 we prove, however, that the A-system associated with a rich, rank(iV) > 
frame allows no strictly hyperbolic solutions. 
• rank(iV) = n — 1. In this case the only solution to (|2.15p is x = 0, that is, all A* are equal. 
According to Proposition 12.11 it follows that all eigenvalues are equal to a common constant: 

rank(iV) = n — 1 =>• A-system has only trivial solutions. 

In particular, if the A-system admits a strictly hyperbolic solution, then necessarily rank(iV) < 
n — 1. It is a non-obvious fact that there are cases where all solutions to f|2.12|) - (|2.13[) are 
trivial; explicit examples are provided by Examples 15.21 15.41 15.51 and 15.101 

On the other hand we observe that the condition rank(iV) = n — 1 does not characterize 
the cases where the A-system (12 . 1 2 H - ( |2 . 1 3 has only trivial solutions. In other words, it may 
be that the only solutions to (|2.12p - (|2.13l) are the trivial ones, while rank(iV) < n — 1; for 
a concrete example see Example 15.21 and 15.41 Thus: 

A-system has only trivial solutions 7$- rank(iV) = n — 1. 

We summarize our findings: 

Proposition 2.2. Consider Problem[l\ for a given frame {Ri(u), . . . , R n (u)}. Then the \-system 
may be formulated as an algebraic-differential system $2.12\) - [2.13\) for the eigenvalues A*. Absence 
of the algebraic constraints 12.15]) implies that the frame is rich, but not vice versa. A maximal 
number of n — 1 independent algebraic constraints 112.1 Sty implies that the \-system admits only 
trivial solutions, but not vice versa. 



3. Systems of three equations 

In this section we present a complete breakdown of the possible solutions of Problem [T] for a 
given frame of vector- fields in R 3 . When n = 3 the A-system (|2.12p - (|2.13p consists of six linear 
PDEs and three linear algebraic equations. For concreteness we record these; the PDEs are 



(3.1) 


n(A 2 ) 


= ri^-A 2 ) 


(3.2) 


n(A 3 ) 


= rjji^-A 3 ) 


(3.3) 




= rJ^-A 1 ) 


(3.4) 


r 2 (A 3 ) 


= r! 2 (A 2 -A 3 ) 


(3.5) 




= r i3(A 3 - A 1 ) 


(3.6) 


r 3 (A 2 ) 


= r 2 3 (A 3 -A 2 ) 



while the algebraic constraints may be written as 



(3.7) 



Nx 



(r 2 



r 1 

1 32 



31 



r 3 
1 12 



rf 3 ) 



(r 



-r 1 

1 23 

r 2 

1 13 



3 
21 



r? 2 ) 



0. 



where x 2 
I 



A^ 



A 1 and x 3 = A 3 — A 1 . There are three possibilities depending on rank(iV): 
rank(iV) = 0. There are no algebraic constraints; the eigenvectors are pairwise in involution, 
and any corresponding system of conservation laws (jl.ip is rich. The analysis in Section 
14.21 which applies to systems of any size, demonstrates that the A-system always has many 
non-trivial (in particular, many strictly hyperbolic) solutions in this case. 
II: rank(iV) = 1. In this case (J3TTJ) imposes a single linear relationship among the eigenvalues. 
This case is more involved and there are several possibilities in terms of how many constants 
and functions determine a general solution. The analysis is detailed in Section 13.11 below. 



Ill: rank(iV) = 2. In this case there are only trivial solutions A 1 

Section [S] provides examples for each type of behavior. 
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A" = A J 



constant. 



3.1. Case II: a single algebraic relation. Using (| 1 . 15[) we rewrite the algebraic relations f|3. T|) : 

(3.8) 4A 1 = Tl^-T^X 3 , 

(3.9) c 2 31 x 2 = rl^-rfgA 3 , 

(3.10) c 3 21 x 3 = rf^A 1 - rf 2 A 2 . 

By assumption the rank of the system (|3.7p is 1, whence the three equations are all equivalent to 
a single non-trivial algebraic condition (unique up to non-vanishing scalings) 

(3.11) aiA 1 +a 2 A 2 + a 3 A 3 = 0, 
where necessarily 

(3.12) a 3 = -(a 1 + a 2 ). 
There are therefore two sub-cases to consider: 

• Ila: all three A 1 appear in (|3.1ip with non-zero coefficients, 

• lib: only two of three A* are involved in (|3. 1 1|) with non-zero coefficients. 

In either case it may be that the only solutions are trivial. To analyze non-trivial solutions we 
employ the Frobenius integrability theorem in case Ila, while case lib requires the more general 
Cartan-Kahler integrability theorem. It will turn out that the number of non-trivial solutions differ 
in the two situations. 

3.2. Subcase Ila: All three A 1 appear in the unique algebraic relation. We first recall the 
relevant formulation of the Frobenius integrability theorem. 

Definition 2. Let U C M. m , V C M. n be open sets, and assume that : U X V — > R, 1 < j < m, 
1 < i < n, are smooth functions and {Yi, . . . ,Y m } is a frame on U. Then the first order PDE 
system 

(3.13) Yj(v % ) = g Z j(x, v(x)) for the unknown n-vector v(x) = (u 1 (x), . . . , v n (x)), 
is called a Frobenius system (in n unknowns onU x V). 

That is, a Frobenius system prescribes all first derivatives of all the unknowns, see 



Theorem 3.1 (Frobenius Integrability Theorem - Frame Version). Suppose $3.13]) is a Frobenius 
system in n unknowns onlA x V that satisfies the following integrability conditions as identities in 
(x,v) eWxV: 



in 



(3.14) Y k (<?;) - Yj (gl) = <Aj9i for 1 < j < k < m, 1 < i < n . 

i=i 

Here the structure coefficients a l k -, 1 < j, k, I < m, are given by [Yk,Yj] = J2h=i a kjYi> an d ^k{g l j) 
denote the total derivatives obtained by using the equations in $3.13}) : Yk{g % j) := Yj c (g i j)+ yJu9j) '9k ■ 
Then, for a fixed point (x,v) G U x V (and under suitable smoothness conditions on g 1 -, Y^), the 
system 13.13\) has a unique local solution v(x) defined for x near x. Furthermore, these solutions 
foliate a neighborhood of (x, v) as v varies over V, and the general solution to f3.13\) depends on n 
constants. 

In the case Ila it turns out that the A-system can be rewritten as a Frobenius system and we have: 

Theorem 3.2. Assume n = 3 and that the \-system contains a single algebraic constraint \3.11\i 
frank(7V) = 1). Consider the case Ila where all three X 1 appear with non-vanishing coefficients in 
13.11\) . Then, after elimination of one of the unknowns, the PDEs in the X-system, reduce to a 
Frobenius system for two unknown functions of three variables. If the corresponding compatibility 
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conditions \3. 14\ ) are satisfied as identities then the general solution to the X-system \3. l\) - [3Td\) . 
\3. l\ j depends on two constants. Otherwise, the only solutions are trivial solutions. Finally, both 
situations can occur. 

Proof. We return to (|3.1ip . where we assume that all three X 1 are non- vanishing. By relabeling 
indices if necessary we may assume that (|3. 11 [) is proportional to (|3.8p . and hence c\ 2 7^ 0, T\ 2 ^ 
and / 0. We solve <^B) for A 1 : 



(3.15) 



1 \3> 



32 



and use this to eliminate Ai in the differential equations (I3.1|) - (l3.6p . By using repeatedly that 
c ij = ^ij ~ ^jv sucn tnat r m(-pr) = r m {-Jr) for all i,j,k,m, we obtain the PDEs 



(3.16) 

(3.17) 

(3.18) 
(3.19) 
(3.20) 
(3.21) 



r 2 ('EpA (A 2 _ A 3 } + Ep r2 (A 2 ) 

c 32 / c 32 



r 3 (A 2 -A 3 ) + ^r 3 (A 2 ) 

c 32 / c 32 



n(A 2 ) 
n(A 3 ) 



%r 2 (A 3 ) 

c 32 



r 2 r 1 

1 21 1 23 / \2 



(A - A 3 ) 



32 



r 3 r 1 

i3li32 (A 2_ A 3 )5 



32 



r 1 r 1 

il2£23 (A 3 _ A 2 } 



32 



r 2 (A 3 



r 1 



lfr 3 (A a 



32 



r 3 (A 2 



r 3 2 (A 2 -A 3 ), 

L_13hu (A 3 _ A 2 } 

C32 

r 2 3 (A 3 -A 2 ), 



Since r 32 , r 23 / 0, we can solve (^T8D and (1330]) for r 2 (A 2 ) and r 3 (A 3 ) by using (|3.19p and (|3.2ip : 



(A 2 ) = Sfk (A 2 _ A 3 } ) ri(A 3 } = £|ll32 (A 2 _ A 3 } 



32 



32 



(3.22) r 2 (A 2 ) 



r 23/ r 3 r l \ C 32 f^32 
^Ti L 32 - 1 12 J - fa" r 2 I 



32 
1 32 /-pi 

■^r\ L 13 

1 23 



-32 



1 23 



32 
1 23 



32 



(A 2 - A 3 ) , 



(A 2 - A 3 ) , 



r 2 (A 3 )=r 3 2 (A 2 -A 3 ), 
r 3 (A 2 ) = -r 2 3 (A 2 -A 3 ). 



This system specifies the derivatives of the two unknown functions A 2 and A 3 along all three vector 
fields ri, r 2 and r$. Hence the system is of Frobenius type. For simplicity we write the system as 

(3.23) n{\ s ) = (j)t(u)(\ 2 - A 3 ) for i = 1,2, 3 and s = 2, 3, 

where (fij are known functions of T's, given by the right-hand sides in (|3.22p . According to the 
Frobenius Integrability Theorem this system is integrable provided 

3 

(3.24) [n ,rj]=Y^ 4 for 1 < i < j < 3 , 

fe=l 

where Cy = T^- — and the left-hand side is computed by using the equations (|3.23|) . A calculation 
reduces (pT2lj) to: 
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(3.25) 



r i (^)-r i (0|) + 



where 1 < i < j < 3 and s = 2, 3. 
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) (A 2 - A ) 



,fe=i 



4^ 



(A 2 - A 3 ) , 



These conditions are satisfied if A 2 = A 3 , in which case the system (|3.22p implies that A 2 = A 3 is 
a constant. Equation (|3.15p then shows that A 1 = A 2 , and Proposition 12.11 implies that the solution 



in this case is trivial: A 



A 2 = A d 



constant. 



For a non-trivial solution to exist the following six conditions must hold: 



(3.26) 



(3.27) 



r; 



fc=i 

m 

k=l 



1 < i < j < 3, 



1 < i < j < 3. 



Example 1 5 . 1 1 and Example 15 . 21 show that these compatibility conditions may or may not be satisfied: 
they must be checked for each case individually. If the compatibility conditions are met then, 
according to the Frobenius Theorem, the general solution to the A-system depends on two constants. 

□ 



3.3. Subcase lib: Exactly two A* appear in the unique algebraic relation. This case is 
more involved than Ha: the A-system does not reduce to a Frobenius system and the non-trivial 
solutions must be analyzed by using the more general Cartan-Kahler theorem. This requires a 
reformulation of the A-system as an Exterior Differential System (EDS). Using the terminology and 
notation of [7] and [20] we give the following formulation of Cartan-Kahler theorem (see Theorem 
7.3.3 in [20] and discussion on p. 87 in [7]). In order to apply this theorem we need to assume that 
the given frame is analytic. 

Theorem 3.3. (Cartan-Kahler Integrability Theorem). Let E^, < k < n, be a flag of integral 
elements at a point p for an analytic EDS on an (n+s)- dimensional manifold, with dimi^ = k, and 
such that Ek is Kahler regular. Then there exists a smooth n- dimensional integral manifold S whose 
tangent space at p is E n . Furthermore, let H(E^) be the polar spaces of E^ and Cf- = codimff(£ , / c ) 
for < k < n— 1, with c n = codim.E n . Let x = (x , . . . , x n ) and y = (y , . . . , y s ) be local coordinates 
around p, chosen so that E^ is spanned by gfr> • • • ) gfif ; E n is annihilated by dy l , . . . ,dy n , and 
H{Ek) is annihilated by dy 1 , . . . ,dy Ck , for < k < n — 1. Then S is defined in the coordinates 
(x, y) by analytic equations of the form y a = F a (x 1 , . . . , x n ), a = 1, . . . , s. 

More precisely, letting x k = x k (p), then in a neighborhood of p the integral manifold S is uniquely 
determined by the following initial data: 



(3.28) 



f° 



F a (x l ,x 2 ,... 
F^x^x 2 ,. . . 
F a (x 1 , x 2 , x 3 , 



,x' 



for < a < Co , 
for c < a < ci, 
for c\ < a < C2 , 



f a (x\x 2 ,...,x n ) := F a {x 1 ,x 2 ,...,x n ) 



for c„_i <a<c n . 



Namely, for each k = 0, . . . , n and a such that c&_i < a < c^, let f a (x 1 , . . . , x k ) be an arbitrary 
analytic function of k variables, such that \ f a {x 1 , . . . ,x k ) — y a (p)\ is sufficiently small in a neigh- 
borhood of (x 1 , . . . ,x k ) (where we define c_i := and by a function of zero variables we mean a 
constant). Then there exists a unique analytic integral manifold y a = F 01 ^ 1 , . . . , x n ), a = 1, . . . , s 
that satisfies the initial conditions \3. 28\) . 

We now return to the A-system. To structure the presentation we first relabel indices (if neces- 
sary) and make the assumption 

(A) the unique algebraic relation (13.11h does not involve A 1 , i.e. a\ = 0. 
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Theorem 3.4. Assume n = 3 and that the X-system contains a single algebraic constraint \3.ll\i 
(rank(iV) = 1). Consider the case lib where exactly two X 1 appear with non-vanishing coefficients 
in \3.11\) . and assume without loss of generality a labeling such that the assumption (A) holds. 

IfT\i ^ r 21 then the only solutions to the X-system are trivial solutions A 1 = A 2 = A 3 = constant, 
while if r 31 = then the general solution depends on one arbitrary function of one variable and 
one arbitrary constant. Finally, both situations can occur. 

Proof. The algebraic relation (|3,lip is equivalent to (|3.8|) - (j3.10[) . and from (A) it follows that 

(3.29) 4 = 1^-1^3 = 0, r| 1 = 0, and 11 = 0. 
Thus the algebraic relations (|3,8p - (|3.10p reduce to: 

r 32 (A 2 -A 3 ) = 
r 2 3 (A 2 -A 3 ) = 

T 3 2 (A 2 -A 3 ) = 0. 

By assumption this system has rank 1 such that the algebraic equations in the A-system are satisfied 
if and only if A 3 = A 2 . Using this relation to eliminate A 3 in the PDE system (I3.ip - (|3.6p we obtain: 

(3.30) n(A 2 ) = ri^-A 2 ) 

(3.31) n(A 2 ) = rf^-A 2 ) 

(3.32) r 2 {X l ) = ri 2 (A 2 -A 1 ) 

(3.33) r 2 (A 2 ) = 

(3.34) r 3 (A 1 ) = r} 3 (A 2 — A 1 ) 

(3.35) r 3 (A 2 ) = 0. 

It follows immediately from the first two equations that if r 31 7^ r 21 , then A 1 = A 2 = A 3 , and thus, 
by Proposition 12.11 the only solutions are trivial: A 1 = A 2 = A 3 = constant. On the other hand, if 
A 3 = A 2 , such that the algebraic equations in the A-system are satisfied, and in addition 

(3.36) I* = ri , 

then the A-system reduces to the sub-system f|3.31j) - (|3.35[) of five PDEs on M 3 for two unknowns 
A 1 and A 2 . As before, A 1 = A 2 (= A 3 = constant) provides a trivial solution to the A-system by 
Proposition 12.11 Example 15.41 shows that such A-systems exist. 

To analyze non-trivial solutions we now assume, in addition to (|3T29j) and (|336|h that A 1 / A 2 = 
A 3 . We note that the sub-system (|3.3ip - (|3.35p is not of Frobenius type, and instead we apply the 
Cartan-Kahler theorem. This requires some preliminary calculations. 

We start by verifying the integrability conditions corresponding to equality of all mixed 2nd 
derivatives that can be computed from (I3.3ip - (|3.35p : [r 3 ,r2](A 1 ) = l^c^r^A 1 ), [?"2, ri](A 2 ) = 
C2i r fc(A 2 ), [r 3 ,ri](A 2 ) = ^ c 3 i r fc(A 2 )> and [?" 3 , r2](A 2 ) = ^ c 3 2 r fc(A 2 )> where the summations are 
over k = 1, 2, 3. Taking into account the assumptions A 1 7^ A 2 and c\% = 0, these reduce to: 

(3-37) r 3 (ri 2 )-r 2 (r} 3 ) = + c 3 ^ 

(3.38) r 2 (T 2 21 ) = T 2 ^ 

(3-39) r 3 (T 2 21 ) = T 2 ^. 

Lemma 3.5. Due to symmetry U.15\) and flatness i!.16\) . and the assumptions h3. 29\) and h3.36]) . 
the compatibility conditions (f3. 31 ), \3. 38\) . \3. 39\) are all satisfied as identities. 

Proof of Lemma \3.b\ We only include the most involved case of (|3.37p . Equation (|1.16p with m = 3, 
k = 1, i = 2, and j = 1, after simplifications due to (|1.15p . reads: 

r 3(ri2) — r i(r 3 2) = (riir^ — r^ 3 rj 2 ) + (rj 2 r| 2 — r 3 2r 2 2 — c 2 3 r 22 ) + (r^ 3 r| 2 — r 33 r?2 — c ? 3 r 3 2) • 
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Similarly, equation (|1.16p with m = 2, k = 1, i = 3, j = 1 gives: 

r 2(ri 3 ) — n.(r 23 ) = (r^r^ — r^r^) + (r| 2 r 23 — r 22 r 2 3 — c 2 2 r 23 ) + (r|[ 3 r 23 — r^r^ — ci 2 r 33 ) . 

Subtracting these and keeping in mind (I1.15|) and assumptions (|3.29H and (13.36p . we get 

r 3(r 12 ) — r 2(r 13 ) = (r^r^ - r 23 r 2 2 — c 2 3 r 22 + rj 3 r| 2 — r^r^ — ci 3 r 23 ) 

/pi -p2 -pi -p2 2 -pi _i_ pi p3 pi p3 3 pi \ 

~~ \ L 12 1 23 — 1 22 1 13 ~~ C 12 L 23 + 1 13 1 23 — 1 23 1 13 ~~ c 12 i 33) 

— r 1 ^ 2 r 1 r 2 _i_ r 1 r 2 _i_ r 1 ,-> 3 r 1 r 3 _l r 1 r 3 

— 1 12 c 32 23 21 ~~ ^ 22 1 31 "+" 1 13 c 32 ~~ 1 33 1 21 ' 1 23 1 31 

= rj 2 c| 2 + rj 3 c| 2 . 

which is (|3.37p . Similar calculations verify (|3.38p and (|3.39p . □ 

Returning to the proof of Proposition 13.41 we proceed to analyze the EDS I associated to the 
PDE system (I3.3ip -( [3.35p . which is differentially generated by the 1-forms 

e 1 = dA 1 -s£ 1 -r} 2 (A 2 -A 1 )£ 2 -ri 3 (A 2 -A 1 )^ 
e 2 = d\ 2 - r^^A 1 - a 2 ) i 1 , 

on M := M 6 with coordinates u , u 2 , u 3 , A , A 2 , and s, where s represents ri(A 1 ). Lemma 13.51 
amounts to the fact that the essential torsion of I vanishes identically. A direct computation 
(making use of Lemma 13. 5| c 32 = and T 31 = T 21 ) shows that 

dO l = irAi 1 mod{#\6> 2 } 
d6 2 = mod^ 1 ,^ 2 }, 

where 

vr = -^+[-5(c} 2 + rl 2 ) + (A 2 -A 1 )(r 1 (rl 2 )-rl 2 r 2 2 -rl 3 c? 2 )]£ 2 

(3-40) + [ - s(c\ 3 + r} 3 ) + (A 2 - A 1 )^! (T\ 3 ) - T\ 3 T 3 3 - T\ 2 c 2 3 )} £ 3 

=: -ds + A 2 £ 2 + A 3 £ 3 . 

Thus, the 1-forms £ 1 ,£ 2 ,£ 3 ,dX 1 ,d\ 2 ,n provide a coframe on M 6 , and the EDS I is algebraically 
generated by 6* 1 , 6> 2 , and ir A £ . 

The remaining parts of the proof consist in: (1) describing the variety of 3-dimensional integral 
elements of 2 that satisfy an independence condition, (2) choosing a flag of integral elements, and 
computing corresponding polar spaces, and (3) applying the Cartan test and the Cartan-Kahler 
theorem to determine the set of solutions. 

(1) The variety of integral elements: let E 3 6 G 3 (TM| P ) be an integral element of I at the 
arbitrary point p £ M = M 6 , where G 3 (TM| p ) denotes the Grassmannian manifold of 3-dimensional 
subspaces of the tangent space T p M. The 3-plane E 3 is required to satisfy the independence 
condition 

(3.41) £ l A£ 2 A£ 3 \e s 
We let {ei, e 2 , e 3 } be a basis of E 3 such that 

(3.42) ^(ej) = Sij i,j = l, 2, 3. 

The following are then necessary and sufficient conditions for E 3 to be an integral element of Z: 

1 (e i ) = O, d 2 ( ei ) = and tt A ^(e*, ej) = , where i,j = 1, 2, 3 . 
Due to conditions (I3.42p these are equivalent to 

(3.43) 9 1 (e i ) = 0, 9 2 ( ei ) = 0, for i,j = 1,2,3 and vr(e 2 )=0, 7r(e 3 ) = 0. 
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Thus the variety V^{T p ) C G 3 (TM\ p ) of 3-dimensional integral elements is defined by 8 independent 
linear equations (|3.43[) . We conclude that 



codim V 3 (T| P ) = 8, and so dimV 3 (J| p ) = dimG 3 (TM| p ) - codim V 3 (J| P ) = 9-8=1. 
Explicitly we can write 

^ =n+p\£r+p 2 i ^i+a i -§- s , for i = 1,2,3, 

where p\ , pf , aj , i = 1, 2, 3, are coordinate functions on an open subset of G 3 (TM\ p ) on which the 
independence condition (|3.41|) is satisfied. The conditions (|3.43p imply that V 3 (T) is defined by the 
equations 

(3.44) p\ = s, p^Tl^-X 1 ), p\ = rj; 3 (A 2 - A 1 ) , p\ = ^(A 1 - A 2 ) , 
and 

(3.45) P 2 2 = 0, P 2 3 = 0, a 2 = A 2 , a 3 = A 3 , 

where A 2 and A 3 are defined in (|3.4(jp . The variety V 3 (Z) is therefore parametrized by a unique 
coordinate function a\. Equations ()3.44p are of constant rank at every point p G M, and therefore 
all integral elements are ordinary at every point of M. 

(2) Flag of integral elements: By setting the arbitrary parameter a\ = we specify a particular 
integral plane E 3 = spanjei, e 2 , e 3 }, where 

ei = n + ^ + rl^A 1 -^)^, 

e 2 = r2 + rl 2 (A 2 -A 1 ) 5 f r + A 2 £, 
e 3 = r 3 + rl 3 (A 2 -A 1 )gf r + A 3 f . 

We define a flag of integral elements: 

(3.46) E = {0} C Ei = span{ei} C E 2 = span{ei,e 2 } C E 3 = span{ei, e 2 , e 3 } , 

A calculation, using that 7r(ei) = 0, ^(ei) = 1, and ^ x (e 2 ) = 0, shows that the corresponding polar 
spaces are 

H(Eo) = {6\e 2 }^, 

H(E 1 ) = {v G H(Eq)\tt A l l (e\,v) = 0} = {9 1 , 9 2 , ir} ± 

H{E 2 ) = {v G H(E )\tt Ai 1 (e 1 ,v) = and 7r A ^ 1 (e 2 , f ) = 0} = {9 1 ,6 2 ,7r} ± . 
Therefore 

Co := codim H(E ) = 2 , ci := codim H{E{) = 3 , c 2 := codim H(E 2 ) = 3 . 

(3) Applying Cartan test and the Cartan-Kahler theorem: Since Co + c\ + c 2 = codim V 3 (X) = 8 the 
equality in the Cartan's Test for involutivity holds. Therefore each E^ in the flag (|3.46p is Kahler 
regular (see Theorem 7.4.1 |20|). 

According to the Cartan-Kahler Theorem 13.31 the general solution depends on 2 constants that 
prescribe the values of A 1 and A 2 = A 3 at an initial point u and one arbitrary function of one variable 
that prescribes the directional derivative s = ^(A 1 ) along a curve. This arbitrary function absorbs 
the arbitrary constant that prescribe the value of A 1 at u and thus effectively the general solution 
depends on one arbitrary function of one variable and one arbitrary constant. The A-system in 
Examples 15.31 has solutions of this type. This concludes the proof of Proposition 13.41 □ 
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4. Rich systems 



In this section we consider Problem Q] for rich systems where we are given a coordinate system u 
and a frame {R\(u), . . . , R n (u)} satisfying Definition [TJ This corresponds to searching for systems 
(jl.ip whose eigencurves are the coordinate curves of some system of coordinates w l (u), . . . ,w n (u) 
(the Riemann invariants) on ]R n . This leads to a slight reformulation of Problem [IJ and a theorem 
of Darboux provides a complete solution for any dimension n. 

The reason why general rich systems can be completely analyzed in this way is the fact that 
all the algebraic constraints in this case either are absent, or always impose equality of pairs of 
eigenvalues. This is in contrast to general (non-rich) systems for which the algebraic conditions 
may be more complicated; see Section 13.11 

4.1. The A-system in Riemann invariants. We recall that richness may be formulated in several 
equivalent ways. (For the setting where a system (jl.ip is given, see Chapter 12 in [35] or Sections 
7.3-7.4 in [E].) m particular, it follows from Definition Q] and the Frobenius theorem that the frame 
{R±(u), . . . , R n (u)} is rich if and only if there is a change of coordinates 

u p(u) = (w l (u), . . . ,w n (u)) with Vw l (u) ■ Rj(u) 

The ^-coordinates are referred to as associated Riemann invariants. These are not unique and we 
assume that we have fixed one choice of the map w = p(u). In this case we may scale the given 
vector fields Ri(u) according to the normalization 

(4.1) Vw\u)- Rj{u) = 5), 

which we assume throughout this section. A calculation (see [12| Section 7.3) then shows that 
these scalings render {Ri(u), . . . , R n (u)} a commutative frame: all structure coefficients c l - k in (jl.5p 
vanish. Furthermore, according to the normalization (|4.1j) . the corresponding left eigenvectors 
L l (u) in (|1.17p 9 are given by L l (u) = Vw l (u), such that the matrix L(u) is a Jacobian with respect 
to u: Ldu = dw. 

In ^/-coordinates the A-system (|2.12p - (I2.13P becomes 

(4.2) = Z^-kP) for l<i#j<n, = JL) 

(4.3) Zt(K j - k { ) = for 1 < k ^ i < j ^ k < n, 
where 

(4.4) «*(«/) := A* o p-\ w ) and Z^(w) := I* o p" 1 ^) . 

Symmetry (|1.15p and flatness (|1.16p of the connection V imply the following properties of the 
Christoffel symbols Z^(w): 

(4.5) Zg = Z% (symmetry), 

n 

(4.6) d m (4) - d k {Zl m ) = (ZiA - Z\ m Z\ k ) , Vt, j, k, m (flatness). 

t=i 

Problem Q] in the rich case thus takes the following form: 

Problem 2. (Rich frame) With the same notation as in Section \1.3l assume that the given frame 
{R\(u), . . . , R n (u)} is rich, and let w = (w 1 ,...,w n ) be associated Riemann invariants. Assume 
the normalization |^.i[ ) ; and define the connection coefficients Z^- by Then, determine the 

set of solutions n l {w), . . . , K n (w) of the X-system in Riemann invariants KJJty - fiJlfy . 
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/ =0 ifi^j, 
I/O if i = j . 



Below we give a complete answer to this problem by applying a theorem of Darboux. For 
completeness we include a precise statement of this result. 

Theorem 4.1. (Darboux [13] ) Given a system of first order PDEs with dependent variables v = 
(v 1 , . . . , v m ) and independent variables x = (x l , . . . , x n ). Assume that each equation is of the form 

(4.7) ^ = f ih (v,x), 

where the given maps v \— > fih(v, x) are C 1 -smooth, uniformly in x in a neighborhood of a given point 
xq = (xq, . . . ,Xq). Assume that the system prescribes compatible second order mixed derivatives in 
the following sense: 

(C) Whenever v l is such that two different first derivatives and are given by the system 
d^. 7j ), the equation 

d d 

g^{fih{V,X)) = —(fafax)) 

(after taking the total derivative on each side), contains only first order derivatives which 
are prescribed by {4- 7\ ), and substitution of {4- 7j ) for the first derivatives that appear gives 
an identity in x and v. 

Suppose a dependent variable v l appears differentiated in \4- 7[ ) with respect to x 1 , . . . , x hp . Then, 
letting x, denote the remaining independent variables, we prescribe a real-valued C 1 -smooth function 
<p l (x) and require that 



(4.8) v\x x ,...,x n ) 



X h l=X n 1 ,...,X h P=X n P 



We make such an assignment of data for each dependent variable that appears differentiated in 
\4- ?| ). Then, under the compatibility condition (C), the problem \4- 7| ) - J^.<§| ) has a unique local 
C 2 -solution. 

Remark 4.1. The initial data ( f-^.ffp provide the values of each dependent variable in the directions 
where its derivatives are not prescribed by the system {4- 7\ ). In particular, if all partial derivatives 
of a variable v 3 are prescribed by 7| ), then ip 3 (x) is taken to be a constant. 

Our formulation of Theorem \4-l\ covers several situations of both over- and under- determined 
systems. In [13] the different cases are stated in three separate theorems, one of which is the PDE 
version of the Frobenius theorem (all first partials of all dependent variables are prescribed) . The 
theorem is proved by a suitable Picard iteration and requires only C 1 -regularity. 

4.2. Rich systems with no algebraic constraints. We first consider the situation when Zfj = 
for all triples of distinct i, j, k, such that the algebraic constraints (14. 3D are trivially satisfied, i.e., 
rank(iV) = 0. We verify that the compatibility conditions (C) in Theorem 14.11 are satisfied in this 
case, and that the general solution of the system (14. 2p depends on n arbitrary functions of one 
variable. In particular, there are many strictly hyperbolic solutions of Problem [2] in this case. 

We also show that if the given vector fields Ri (in addition to forming a rich frame) are orthogonal, 
then the corresponding A-system necessarily belongs to this unconstrained case. Two concrete 
examples of rich and algebraically unconstrained systems are the class of 2 x 2-systems, and the 
class of n x n-systems with constant eigenfields, see Examples I5.6[ I5.7[ and 15.81 

Theorem 4.2. Given a C 1 -smooth, rich frame {Ri(u), . . . ,R n (u)} in a neighborhood ofw £ M n . 
Let p{u) = (w 1 (u) , . . . , w n (u)) be associated Riemann invariants and assume the normalization 
J^.i[ ). Let the connection coefficients Z\- be defined by h4-4\ ) ind assume that Zf- = whenever 
i 3 ^ k ^ i. Then, for given functions ipi, i = 1, . . . ,n, of one variable, there is a unique local 
solution k 1 (w), . . . , K n (w) to the X-system \4->ty with 

K^W 1 , . . . ,W i -\w i ,W i+ \ ...W n ) = tp^W*). 
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Proof. The A-system (|4.2[) consists of n(n — 1) PDEs that prescribe, for each j = 1, . . . , n, all first 
partials of n 3 except djK 3 . The conclusion of the theorem thus follows from Darboux's Theorem 14. II 
provided the system f|4.2j) satisfies the compatibility conditions (C). That is, for each j = 1, . . . , n 
the equalities 

dkd m K j = d m d k K j 

should hold as identities for all k ^ j, m ^ j and m ^ k, when the first derivatives of k's are 
substituted for from (|4.2p . This leads to the compatibility conditions 

(4.9) (d m z] k -d k z] m )^ + {z] m zz k + z] k zi m -z] m z] k -d m z]^K k 

~ \ Z jm Z mk + Z jk Z km ~ Z jk Z jm ~ ^k Z jmj K m = , 

for all distinct j, m, k. We verify these conditions by showing that the coefficients of k 3 , K m and 
K k vanish identically due to (|4.5|) and the flatness condition (|4.6|) . We first substitute i = j in (|4.6|) 
to obtain 

n 

9m Z j k ~ d k Z 3 jm = ^ { Z tk Z jm ~ Z tm Z jk) > 
i=l 

which vanishes since k, m, j are distinct and since, by assumption, = for distinct i\, 12, 13- 
This shows that the coefficient of k 3 in (j4.9j) is zero. The arguments for K fc and K m are similar and 
we only consider the coefficient of K k . Interchanging k and i in (|4.6p and setting z = j yields 

n, 

5 m Zj fc - djZ J mk = ^2 (Z J jt Z^ k - Z J mt Zj k ) . 
t=i 

Again, using (|4.5p . that Zl\ = for distinct i±, 12, and that k, m, j are distinct, we obtain 

OW^jfc — V ml: + ^jk^mk - ^mj^jk ■ 

□ 

4.2.1. Rich, orthogonal frames. Assume now that the given rich frame {R\(u), . . . , R n (u)} has 
the additional property that it is orthogonal: Ri-Rj = if i 7^ j. That is, we search for systems (ll.lj) 
whose eigencurves are the coordinate curves of an orthogonal coordinate system (w; 1 (w), . . . ,w n (u)) 
on M n . In this case we show that the connection components Z\- necessarily vanish whenever i, j, k 
are distinct, such that Theorem 14.21 applies. 

We define a matrix S(w), whose components SUw) := Rj o p~ 1 (w) are the pull-backs of compo- 
nents of the matrix R under p -1 . The connections one- forms (see (|1.8p and (|2.10p ) are then given 
by: 

n 

(4.10) /4 ■- J2 Z L du}k = ( s ~ ldS ) j i ■ 

k=l 

By orthogonality we have that 

(4.11) 5 T 5 = diag[|5i| 2 ...|5 n | 2 ] such that S' 1 = diag [(Sil" 2 . . . \S n \~ 2 ] S T . 
It follows that the connection matrix p is given by 

H = S^dS = diag [|5i|" 2 . . . \S n \~ 2 ] S T dS. 
Differentiation of (|4.1ip gives 

d(S T )S + S T dS = diag[d|5i| 2 . . . d\S n \ 2 } , 
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T 



such that 

d|^i| 2 d\S n \ 2 

|Si| 2 ■■■ \S n \ 2 

Component-wise we thus have 



diag 



diagH^I 2 



\S n \ 2 ] fi diag 



(4.12) 



I C.|2 . 



~ J \Si\ 2 
From (J4"7L0|) and P~T2j) it now follows that 



and 



Mi 



2|S, : | 2 



d\Si 



Z h 



I'M 2 

\SA 2 ^ 1 



By symmetry in the lower indices (I4.5h we get, for indices i ^ j ^ I ^ i, that 

2 



(4.13) 

while, at the same time, 
(4.14) 



I'M 2 

I5,;l 2 ^ 



I'M 2 
|S,;I 2 



1 5- 1 2 Jl ' 



Z h 



Z h 



\3l 7 i 
\s-\ 2 3 



\Sil 7 i 



|<SVi| 



It follows from (|4.13p and (|4.14p that = whenever all three indices are distinct, such that 
Theorem 14.21 applies. We conclude that for rich, orthogonal frames the solution of Problem [2] 
depends on n functions of one variable. 



4.3. Rich systems with algebraic constraints. We next consider the more involved situation 
of rich systems with non-trivial algebraic constraints (|4,3p such that rank(iV) > 0: there exist 
distinct i, j, k such that Zf- ^ 0. The algebraic relations (|4.3p then impose the equality k 1 = k? , 
such that there are multiplicity conditions on the eigenvalues. (In particular, there are no strictly 
hyperbolic systems in this "rich & rank(iV) > 0" case.) This makes the analysis of the A-system 
more involved. 

The analysis is further complicated by the fact that there may be additional equalities among 
eigenvalues that are imposed by the PDEs (|4.2p . We claim that, after taking all of these relations 
into account, the A-system (I4.2p - (|4.3p reduces to a first order system of PDEs to which the Darboux 
theorem (Theorem 14. ip can be applied. We have: 

Theorem 4.3. Given a C 1 -smooth, rich frame {R\(u), . . . ,R n (u)} in a neighborhood ofw £ W 1 . 
Let (w (u), . . . ,w n (u)) be associated Riemann invariants and assume the normalization fl^. Let 
the connection coefficients Z^- be defined by \4-4\) > an d assume that there exists at least one triple 
of distinct i, j, k with Z\- ^ 0. 

Then the X-system l4-^ )- (T^\ ) imposes multiplicity conditions on the eigenvalues in the following 
sense. There are disjoint subsets A\, . . . , A SQ C {l,...,n} (sq > 1) of cardinality two or more, 
and such that |j£p-|j5|) impose the equality k 1 = k 3 if and only if i, j £ A a for some a £ 
{1, . . . , so}- Let I = Y2a=i \Aa\ < n an d s\ = n — l. By relabeling indices we may assume that 
{l,...,n}\U2Li4* = {l,...,si}. 

Then, given s± functions of one variable tp 1 , . . . ,tp Sl , and so constants h , . . . , h °, there is a 
unique local solution « , . . . , n n of the X-system ^-Sty - fcJ^) such that 

^(w 1 , . . . ,w l ~ 1 , w\ w l+1 . . . w n ) = fi(w l ) if i = 1, . . . ,si , 

n l (w) = h if i G A a for some a = 1, . . . ,sq. 
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Moreover, diK J (w) = whenever i, j G A a for some a = 1, . . . , so- Here s\ is the maximal number 
of simple eigenvalues, while sq is the maximal number of non-simple eigenvalues, in a solution of 
the X-system. 



This result is a direct consequence of the following two lemmas and Theorem 14.11 The first 
lemma details how the index sets A a are defined, and how to write the A-system (I4.2p -( j4~3j) as a 
system of only PDEs. 



Lemma 4.4. With the same assumptions as in Proposition ^ ■ 3\ the X-system \4-^ - can be re- 
written as a first order PDE system without algebraic constraints. More precisely, there are integers 
so, s\, and disjoint subsets Ai, . . . ,A So C {1, ...,n} as described in Proposition \4-3\ such that 
the following holds. The re-written system involves s\ + So unknowns, denoted ^(w), . . . , k Si (w), 
h}(w), . . . , h s °{w), and has the following form: 



(4.15) for j = : 




if 1 < i + j < si 

if i G A a for some a = 1, . . . , Sq, 



W?(i£-h a ) ifi = l,...,si 

(4.16) for a = 1, . . . , s : dih a = { wf(hP - h a ) if i G Ap and 1 < / a < s 

ifieA a . 

Here the coefficients W!j* are defined by := if j G A a and i G" A a , and these are well-defined 
by the properties of the sets A a . 

Proof. We construct the index sets A±, . . . ,A So by first considering the algebraic conditions (|4.3p . 
and then taking into account the differential equations (|4.2p . 

Each of the algebraic relations in (|4.3p is either trivial ("0 = 0") or non-trivial, imposing equality 
of two eigenvalues. More precisely, if 1 < i 7^ j < n are such that there exists k £ with Zf- 7^ 0, 
then necessarily k 1 = k j . By assumption there is at least one such non-trivial algebraic relation. 
Grouping together the indices i of the unknowns k 1 that must be identical according to these 
relations, we obtain a certain number sq > 1 of disjoint index sets A\, . . . , Ag C {1, . . . , n}: two 
distinct ji, 22 belong to the same A a if and only if one of the relations in (|4.3p imposes the equality 
K n = K n . Clearly, the cardinality of each A a is at least two. For a = 1, . . . , sq we introduce the 
unknowns h a by setting 

h a (w) := K j (w) VjeA a , 
which is well-defined by the definition of A a . Also let 

so 

(4.17) A:= (J A a Si :=n-\A\. 

a=l 

If necessary we relabel the unknowns such that 1 denote the unknowns for which the 

algebraic relations (|4.3p impose no multiplicity constraint. 

We next turn to the PDE system (|4.2p . For each pair with i 7^ j, either none, exactly one, 

or both of i and j belong to A. We list the possible cases together with the corresponding form of 
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the PDEs in each case: 



(a) 


1 < i + j < h : 


cW = 


z]^ 




( b ) 


3a : i G A a and 1 < j ' < s± : 


OiK J — 


Zjiih 


— K J ) 


(c) 


3/3 : j G A/3 and 1 < i < s\ : 


dihP = 


z\^ 


-hP) 


(d) 


3*y 7^ 5 : i £ A~f and j € Ag : 


dih s = 


zU h " 


-h 5 ) 


(e) 


3e : i,j&A e : 


djf = 


. 





At this stage we have used all non-trivial algebraic relations imposed by the algebraic part f|4.3[) of 
the A-system. The issue now is that the cases (c) and (d) may impose further algebraic conditions 
since their left-hand sides are independent of the index j: 

• in case (c): for a given pair with j £ Ap and 1 < i < s±, unless all the coefficients 
coincide as k ranges over Ag, we must impose that k % = hP . If so we add i to the index set 
A 3 and replace the \A B \ PDEs dJiP = - hP) in (c) by the single PDE dji 13 = 0. At 
the same time we substitute hP for k 1 in the remaining PDEs in (a)-(c) in which k 1 appears. 

On the other hand, if = Z^ for all k G Ag, then we define Wj 3 := Zj i: and replace 
the same \Ag\ PDEs by the single PDE dihP = Wf («* - hP). 

• in case (d): for a given pair with i G A~ p j G A5, and 7 7^ 5, unless the coefficients Zj^ 
all coincide as A; ranges over Ag, we must impose that h s = K 1 . We then merge the index 
sets A\ and A s and replace the \A S \ PDEs ^/i 5 = Z^W - h 5 ) in (d) by the single PDE 
dih 5 = 0. At the same time we substitute h s for K 1 in the remaining PDEs in (b)-(e) in 
which h? appears. 

On the other hand, if Z^ = Z 3 ^ for all k G Ag, then we define Wf := Z 3 ^ and replace the 
same \Ag\ PDEs by the single PDE d-h 5 = Wf(k> -h s ). 

We continue this process of identifying unknowns n l and h a that must necessarily coincide (enlarging 
and merging index sets), until no further reduction is possible. Since each reduction consists in 
setting unknowns equal to each other, and since setting all unknowns equal gives the trivial solutions 
to the A-system (Proposition 12 . 1| ) , it follows that no contradiction will be obtained in this manner. 
Also, as there is a finite number of unknowns, the process must terminate after finitely many 
reductions. 

At this point we have obtained a certain number sq > 1 of disjoint index sets A±, . . . ,A SQ C 
{1, ... ,71} such that two distinct j±, J2 belong to A a if and only if (|4.2|) and (|4.3p impose the 
equality n n = k 32 (through the reduction process described above). By assumption no further 
algebraic reduction is possible, and it follows that the sets A a have the following properties: 

(4.19) Z 3 l { =Zf 2i =: W? whenever j x , j 2 G A a and i (£ A a , 
and 

(4.20) Z\j = whenever i G A a , j A a and k g {i,j} . 
We also have that 

(4.21) Z$j = whenever 1 < i ^ j < s\ and k g {i,j} . 
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Each of the sets Ap, (5 = 1, . . . , So, is contained in one of th.6 S6ts a = 1, . . . , so, whence the 
latter have cardinality \A a \ > 2. Setting 



A := A a , and s\ := n — \A\ , 



a=l 

it follows that so is the maximal number of non-simple eigenvalues, and si is the maximal number 
of simple eigenvalues, in a solution of the A-system. If necessary we relabel the indices such that 
1 denote the unknowns for which A4.2[ )- fl4.3j) impose no multiplicity constraints, and we 
introduce unknowns 

(4.22) h a (w) := n j (w) for j G A a , a = 1, . . . , s . 

Similar to above, for each pair (i,j) with i ^ j, either none, exactly one, or both of i and j belong 
to A, and again there are five possible cases of the PDEs in (|4.2[) . We thus obtain the same type of 
system as in (I4.18P except with si replaced by si and A a 's replaced by A a 7 s. Using (14. 19ft we obtain 
the system (|4.15j) - ()4.16p . Finally, from the last equation in (|4,16p it follows that h a is independent 
of the coordinates w l for i G A a : diK 3 {w) = for all i, j G A a . □ 

Remark 4.2. We note that if So = 1 and \A\\ = n, i.e. k 1 = ■ ■ ■ = K n = h 1 , then 15\ )- (4-16\ ) 

degenerate into the trivial system = 0, i = 1, . . . , n. In this case the only solutions to the 
original problem are the trivial ones: k}(w) = • • • = K n {w) = k G M. If I = n, then s± = and the 
system \4- lb^ j is empty. In this case all eigenvalues k 1 necessarily has multiplicity greater than one. 

The PDE system ()4.15|) - (|4.16|) is in the form required by Theorem 14. 1| and we proceed to verify 
the compatibility conditions (C) in this case. 



Lemma 4.5. With the same assumptions as in Lemma \4-4\ the compatibility conditions (C) in 
Theorem \4-l\ are satisfied for the rewritten X-system ^4.15\ )- [4.16\ ). 

Proof. Compatibility conditions for \4-15 ): For j = 1, . . . , s±, we need to verify that 

(4.23) dkdmK 3 = d m dkK? for k / j, m / j and m ^ k. 

There are two main cases depending on whether m and k belong to a common A a or not, and we 
first treat the latter. 

• Case 1. $a such that {m, k} C A a . There are several sub-cases according to whether m or k 
belong to A. However, in each case the compatibility condition takes the form: 



^ {d™ Z \k ~ ®k Z \m) ^ + { Zj jm Z mk + Z jk Z km ~ Z jm Z jk ~ ®m Z \ 



Z \m, Z mk + Z \k Z km ~ Z \k Z \m ~ ®k Z \ m ) K m = 0, 



jk 

i 

jm^mk T ^ jk^km ^jk^jm u k Zj jm 

with the provision that 

K m = h a , Z™ k = W k if m G A a for some 1 < a < s , 



and 



K k = h p , Z\ m = Wl if k G Ap for some 1 < < s . 

We verify that (|4.24j) is satisfied by showing that coefficients of k? , K m and K k vanish identically 
due to the flatness condition (|4.6p . We first substitute i = j in (|4.6p to obtain 

n 

d m { z ] k ) - d k {zij = {Z^Zij - Z^Zlj) . 

t=i 
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By assumption (|4.20p . for distinct k, m, j with j ^ A, we have 

ZL~ = unless t = m or t = j 
Zfrj = unless t = k or t = j , 

such that 

(4-25) d m (4) - d k {Z\ J = Z'l k (Z™ - Z k kj ) . 

By assumption $a such that {m, k} C A a , whence Z° mk = due to assumption (|4.20p . Thus 

(4-26) d m {zj k )-d k {Zj m )=0, 

such that the coefficient of k? in (|4.24|) vanishes. 

Next consider the coefficient of K k in (|4.24|) . Interchange k and i in (|4.6p and set i := j to get 

(4-27) d m (4) - 5, « fe ) = £ (ZjX, " • 



■^mfe = unless t = m or t = k 



Again, by assumption $a such that {m, A;} C A a , and (|4.20p gives 
(4- 

Z*j k = unless t = j or i = k . 
Thus 

yj _ 7m I yj yk yj yj yj yk _ yj ym , 7 j 7 fc 7 jf 7 j 

Om^j k — ^jm^mk + ^jk^km ~ ^jm^jk ~ ^mk^kj ~ ^jm^mk "+" ^jk^km ~ ^jm^jk ' 

such that the coefficient of M fc in (|4.24|) vanishes. The coefficient of n m is treated similarly. 

• Case 2. 3a such that {m,k} C ^4 Q . In this case the compatibility conditions (I4.23h read 

[ dm (zj k )-d k (zj m )](h a -Kn^o. 

We claim that the first factor vanishes, i.e., (I4.26h is valid in this case as well. Indeed, the derivation 
of (|4.25p is valid in this case as well. By assumption 3a such that {m,k} C A a , such that 
Z™j = Zj°. = Z" due to assumption (|4.19p . Thus the compatibility conditions for (|4.15p are met. 

Compatibility conditions for l[4-16\ ): For a = 1, . . . , sq we need to verify that 

(4.29) d k d m h a = d m d k h a for m^k. 

For a fixed a there are several cases depending on whether both, exactly one, or neither of m and 
k belong to A a . 

• Case 1. If both m and k belong to A a then (|4.29|) amounts to the trivial condition = 0. 

• Case 2. Without loss of generality we assume that m G A a ^ k. The compatibility condition 
(|4.29|) then reduces to 

(4.30) = d m d k h a = d m {W£(K k - h a )) = [d m {W£) - W%Z k m ] {K k - h a ) . 

Since | A a | > 1 there is a j / m £ A a , such that W k = Z 3 - k . We will verify (|4.30p by showing that 
the first factor on the right-hand side vanishes: 

(4-31) d m (W£) - W%Z k m = d m (Z0 - Z] k Z k m = . 

Observing that (|4.20p again yields (|4.28[> . we see that (j4.27|) gives 
d m (Zj k ) = Z 3 - k Z km + Z 3 - m Z™ k - Z 3 - m Z 3 - k 

_ yj yk , yj ( ym yj \ _ yj yk 
- A jk Zj km + A jm\ A mk ~ A jk> ~ ^jk^km ) 
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where we have used that Z™ k = Z 3 - k = W k , since j, m G A a . This verifies (|4.31|) . 

• Case 3. If neither m nor k belong to A a then the compatibility condition (|4.29|) reduces to 

(4.32) d m [W?(K k - h a )} = d k [W^(K m - h a )] , 

where W k = Z J - k and = Zj m , for all j G A a , and with the provision that 

K k = hP if keAp(P^ a) and K m = K 1 if m G A 7 (7 / a). 

There are three sub-cases depending on whether both, exactly one, or neither of m nor k belong to 
some Ap with (3 7^ a. To verify the identity (|4.32p we use the same techniques as above: flatness 
(|4.6|) . and (|4.19p - (|4.2ip . The details are similar to the previous cases and are left out. □ 



Combining Lemmas 14.41 - 14.51 with Theorem 14.11 completes the proof of Proposition 14.31 

5. Examples 

5.1. The Euler system for 1-dimensional compressible flow. 

Example 5.1. The Euler system for 1-dimensional compressible flow in Lagrangian variables is: 

(5.1) v t -u x = 

(5.2) u t +p x = 

(5.3) E t + {up) x = 0, 

2 

where v, u, p are the specific volume, velocity, and pressure, respectively, and E = e + \ is the 
total specific energy, where e denotes the specific internal energy. The thermodynamic variables are 
related through Gibbs' relation de = TdS — pdv, where T is temperature and S is specific entropy. 
Using this, the Euler system may be rewritten as 

(5.4) v t -u x = 

(5.5) u t +p x = 

(5.6) S t = 0. 

The system is closed by prescribing a pressure function p = p(v, S) > 0, and we make the standard 
assumption that p v (v, S) < 0. For a given a pressure function the eigenvalues of |5.^[ )-( f5r^) are: 

\ l = -^Ph, A 2 = 0, )? = ^—p- v , 

with corresponding right and left eigenvectors (in (v,u, S) -space and normalized according to Ri ■ 
V ^ 81): ' 

(5.7) i?i = [1, 0] T , R 2 = [- Ps ,0,p v ] T , R 3 = [1, -v 7 ^, 0] T , 
and 



(5.8) L 1 = \ 



1 



1 vs_ 

~Pv ' Pv 



0, 0, 



T 3 _ 1 



1 1 PS. 



The fact that the eigenfields depend only on the pressure allows us to pose the "inverse" problem: 

Problem 3. (Problem^ for the Euler system) For a given pressure function p = p(v, S) > 0, with 
p v < 0, consider the frame given by |5. 7| ), Then: determine the class of conservative systems with 
these as eigenfields by solving the associated X-system i2.12\) - [2.13\) for the eigenvalues A 1 , A 2 , A 3 . 

Remark 5.1. As explained in Remark \l.l\ the two forms \5. l\) - {575]) and (f5.^| )-( f376|) of the Euler 
system are equivalent as far as Problem is concerned. However, they are not equivalent at the 
level of computing weak solutions to Cauchy problems [38] . 
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We shall see that there are two distinct cases depending on whether the frame is rich or not. 
From |5. 7|), 15.8]) . and \2. 8\) we obtain 



21 



r 



2:-! 



o, 



31 



"pl Pvv 



r 1 

1 12 



r 3 

1 32 



2 V p v I v ' 



Thus the PDEs \2.12\) in the X-system are given by 



(5.9) 


n(A 2 ) 


= 


(5.10) 


ri(A 3 ) 


= f^-A 1 ) 


(5.11) 


^(A 1 ) 


= f©,(^-A 2 


(5.12) 


^ 2 (A 3 ) 


= ^f^-A 2 


(5.13) 




= f^-A 3 ) 


(5.14) 


r 3 (A 2 ) 


= 0. 



where r±, T2, r% are the derivations corresponding to R±, R2, -R3. 
algebraic equations \2.13\) of the X-system are: 



r 1 

1 32 



r 3 
1 12 



« (Ps\ 

I \pvlv' 



V 1 

1 23 



r 3 
1 21 



.Pv(PS_\ 
4 V p v J v 



The coefficients appearing in the 
0. 



r 2 

1 31 



r 2 

1 13 



Using these we see that the matrix N in \3. 7j ) is given by 



2 1 

2 1 



By assumption p v < 0, and we thus consider the two cases: 

(a) (^-) v = 0: rank(TV) = and there are no algebraic constraints among the eigenvalues. 
Proposition \2.2\ and Theorem \4-2\ show that the Euler system is rich in this case, and that 
the general solution to Problem^ depends on three functions of one variable. 
(6) (— J ^ 0: rank(iV) = 1 and there is a single algebraic relation among the eigenvalues: 



(5.16) 



A 1 + A 3 



2A 2 



For case (a) the pressure must be of the form p(v,S) = n(£), where £ = v -+- F(S) and F is a 
function of S alone. It follows from 15. 9\) and |5. 14\ ) that X 2 is any function of S alone, and from 
\5.11\) and 115. 12\) it follows that A 1 = A(£,u) and A 3 = B(£,u). The functions A and B should 
then be determined from the remaining equations 15.10\) and A5.13\) : 



A $ - y/-n'(£)A u = a(B -A), B $ + ^-W(£)B U = a(A - B) , 

where a = —j^ 1 = — in^fj- These equations must be solved case-by-case. However, by deriving a 
second order equation for A alone, which in turn will determine B, we see that solving for A and 
B requires two functions of one variable. Together with X 2 (S) these are the three functions that 
determine a general solution in case (a). 

We proceed to explicitly solve the PDE system in the more interesting second case (b). Adding 
115.11]) and K5.12\) . and using 15.16]) . gives r 2 (X 1 + A 3 ) = 0. Applying i5.16\) again yields r2(A 2 ) = 0. 
Together with \5. 9\) and fa5.14\ ) this shows that 

(5.17) A 2 = A (constant). 

To continue it is convenient to express the coordinate frame {d v , d u , ds} in terms of the r-frame: 



(5.18) 



d v = Hn + r 3 ) , 



d u 



= (n - r 3 ) 



ds = ±r 2 + §- v (r l+ r 3 ) 
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A direct calculation, using \5.18\) and \5. 1 6\) . now shows that A := (A 3 — A 1 ) satisfies d u A = and 

d v A=(^)A, d s A=(^)A. 

Integration with respect to v and S, respectively, shows that 

A(v, S) = A 3 - A 1 = D^p~ v 

for a constant D. Invoking \5. 16\) a last time (with A 2 = X), we conclude that the general solution 
of the X- system for non-rich gas dynamics (case (b)) depends on two constants X, C according to: 



A 1 = A - C^p~, 



x- 



X. 



X 3 = X + C^p^ 



5.2. Additional examples with non-rich frames. The case of non-rich gas dynamics in Ex- 
ample 15.11 provides an example of Subcase Ha in Section 13.11 in which the compatibility conditions 
(|3.26[) - (|3.27j) for the Frobenius system (|3.22[) hold identically. The following Example 15.21 of a sys- 
tem in the same category, but with only trivial solutions, shows that these compatibility conditions 
do not necessarily hold as identities and must be checked in each case. This verifies the last claim 
in Proposition! 



Example 5.2. n = 3, NON-RICH SYSTEM with a single algebraic relation of type Ha AND 
ONLY trivial solutions. We prescribe the frame {R±, R 2 , R3} on R 3 by 



The system h3.1 



Ri = [0,0, if 
113. 6\) takes form 



R2 = [0,1,^] 



d 3 X 2 



5, A 3 



^A 1 + u^sX 1 
d 2 X 3 + u^sX 3 
u 3 d 1 X 1 + d 3 X 1 + ^(X 3 -X 1 ) 
u 3 9iA 2 + d 3 A 2 

In this case the matrix N is given by 



R 3 = [u 3 ,0, 1] 

= 0, 

= 0, 

= 0, 

= 0, 

= 0, 

= 0. 



T 





u 3 


u 1 - 








N = 

















and has rank 1. It corresponds to the unique algebraic relation that involves all three eigenvalues: 

(u 3 ) 2 (A 2 - A 1 ) +n 1 (A 3 - A 1 ) = 0. 

The only solution of the above system of differential and algebraic conditions is trivial A 1 = A 2 = 
A 3 = const. 

Next we provide examples of a non-rich 3- frame with a single algebraic relation of type lib. The 
examples show that in accordance with Proposition 13.41 the solution of the A-system in this case 
either depends on one arbitrary function of one variable and one arbitrary constant, or is trivial. 

Example 5.3. n = 3, NON-RICH SYSTEM with a single algebraic relation of type lib 
AND non trivial solutions. Consider a system with eigenvectors: 

Rx = (0,1, Of, R 2 = (1,0, Of, R 3 = (u 2 ,u 3 ,l) T . 

Since [ri,r 3 ] = r 2 , these vector-fields are not pair-wise in involution, and therefore the system is 
not rich. The X-system implies a unique algebraic relations A 2 = A 3 , which involves only two of the 
X's. We are thus in Subcase lib in Section \3. 31 A computation show that = and r|i = are 
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equal and thus according to Proposition \3.4\ we can expect a non-trivial general solution depending 
on one arbitrary function of one variable and one arbitrary constant. Indeed, taking into account 
the algebraic relation the differential equations become: 



n(A 2 ) 
n(\ 2 ) 



diA 1 = 0, 

n^iA 1 + u 3 d 2 A 1 + d 3 X l = 0, 
0, 



d 2 X 2 



diA 2 = 0, 
u 2 dxX 2 



u 3 d 2 \ 2 + d 3 X 2 



0. 



The general solution of this system is A 1 = <p(u 2 — 2u 2 ) , A 2 
constant and ip is an arbitrary function of one variable. 



A, where A is an arbitrary 



Example 5.4. n = 3, NON-RICH SYSTEM with a single algebraic relation of type lib 
and only trivial solutions. Consider a system with eigenvectors: 



R 1 = (u 1 + u 3 ,l,0) T , R 2 = (1,0, 0) T , R 3 = (u 2 ,u 3 ,-u 2 f. 



Since [rx,r 3 ] 



i r i_l_ j^+u^u 3 r2 



r* + 



■ r 3 these vector-fields are not pair-wise in involution, and 



therefore the system is not rich. The X-system implies a unique algebraic relations A 2 = A 3 , which 
involves only two of the A 's. We are thus in Subcase lib in Section Iff. 31 A computation shows 
that = and = 1 are not equal and thus according to Proposition \3.4\ we can expect only 
trivial general solution. Indeed, taking into account the algebraic relation X 2 = A 3 the differential 
equations n(A 2 ) = ^(A 1 - A 2 ) and n(A 3 ) = T^ 1 (A 1 - A 3 ) imply that 



(5.19) 



ri (A 2 ) = (A 2 - A 1 ) and n(A 2 ) = 0, 



which implies A 1 = A 2 . The only non-trivial solution of the A system is therefore A 1 = A^ 
where A is an arbitrary constant. 



A 3 = A, 



We include an example which illustrates the fact that maximal rank of N (rank(iV) = n 
implies that all solutions of the A-system are trivial (see section 12. 3p . 



1) 



Example 5.5. n = 4, system with maximal rank of algebraic constraints (rankA^ = 3). 
The X-system for the following non-rich frame in M. 

R x = (l,0,u 2 ,u 4 ), R 2 = (0,1, u\0), R 3 = (u 3 , 0,1,0), R 4 = (1,0,0,0) 

includes three independent algebraic relations: 

X 3 = X 2 , 

which immediately shows that all solutions must be trivial: A 1 = A 2 
5.3. Additional examples with rich frames. 



A 2 = A 1 , A 3 = A 4 , 



X 3 = X 4 



const. 



Example 5.6. n = 2 SYSTEM. Suppose we want to find the most general system (COP of two 
conservation laws whose eigencurves are hyperbolas and radial straight lines. As a 2 x 2-system it 
is necessarily rich and as Riemann invariants we use w 1 := u 1 u 2 and w 2 := where (u^u 2 ) are 
standard coordinates on M. 2 . Then the matrix of left eigenvectors is 



L(u l ,u 2 



' Vw 1 " 




Vw 2 





TyJ 
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and the matrix of right eigenvectors is R = L 1 . A direct computation shows that Z 2l = ;pr end 
Z\ 2 = 0, such that the X-system {4->ty becomes 

(5.20) d 2 K 1 (w 1 ,w 2 ) = 

(5.21) dm^w^w 2 ) = -L^-k 2 ). 

and A5.21\) then 



By \5. 20\) the first eigenvalue is an arbitrary function of w 1 alone, k 1 
yields 



1 



k 2 ^ 1 , w 2 ) = —= [h(w 2 ) + giw 1 )] , where giw 1 ) 
v w 1 



and h is an arbitrary function of w 2 alone. In accordance with Theorem \4-2\ we see that the general 
solution depends on two functions of one variable. 

As we proved in Section 14.2.11 for a rich and orthogonal frame the algebraic part of the A-system 
becomes trivial. We consider a concrete case in this class. 

Example 5.7. rich orthogonal frame. Suppose we want to find the most general system M.l}) 
of two conservation laws whose eigencurves are the coordinates curves of spherical coordinates (ra- 
dial lines, latitudes, and longitudes) . The corresponding Riemann invariants are polar coordinates 
(r,9,</>). More precisely we let p : R 3 \ {(u l ,0,u 3 ) | u l < 0, z < } -> R + x (-7r,7r) x (0,7r) be the 
change from rectangular Cartesian coordinates (u ,u 2 ,u 2 ) to spherical polar coordinates (r,9,(j)) ; 



(5.22) 

where 



r = \J (u 1 ) 2 + (u 2 ) 2 + (u 3 ) 2 , 9 = arctan 
Defining R by R~ l = Dp we have 

R- 



p{u,v,w) = (r,9,4>), 

( V(n 1 ) 2 + (n 2 ) 2 ' 



and 6 = arctan —r 
v u 1 



u- 



sm V cos < 
sin 9 sin < 
cos 9 



r cos v cos i 
r cos 9 sin < 
—r sin 9 



— r sine/ sin < 
r sin 8 cos $ 




A calculation now yields that Z\ 2 



7 1 
z 13 



7 2 
z 23 



0, z 2 



21 



v 31 



-, Z 3 2 = arctan(#) and, as 



expected, all Z 's with three distinct indices are zero. Thus according to ft4-.2\ )- ffl73\ ), the \-system 
contains only the following set of differential equations: 

dgK 1 = d^K 1 = , 

1 



d r K 



■{K 1 -K 2 ) 



d r K 3 = -{K 1 - K 3 ) 



<v 2 = > 

dgn 3 = arctan 8{k 2 — k 3 ). 



The general solution of this system depends on three arbitrary functions F, G, H each depending on 
a single variable: 

i 



K 



K 



F(r), 
H(6)+g(r) 



+ cos grig | g {r) 



where g(r) := / F(^)d^, 



r 8in( 
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Example 5.8. Constant eigenfields. If the frame {R\(u), . . . , R n (u)} consists of constant 
vector fields R\{u) = Ri then it is rich according to Definition^ Furthermore, all connection 
coefficients Zf- vanish and Theorem \4-2\ applies. As associated Riemann invariants we may choose 
w l (u) = U -u, where each Li is constant. A calculation shows that the general solution of X- system 
{Ij-ty is given by n 1 = ^(w 1 ), . . . , (p n (w n ), where (f> , . . . , (j) n are arbitrary functions of one variable, 
in accordance with Theorem That is, the eigenvalues solving Problem [I] in this case are given 
by \ l {u) = <f) l (L l ■ u), i = 1, . . . , n. 

We have seen an example of rich 3 x 3-system with no algebraic constraints (rankiV = 0) when 
we analyzed eigenvector-fields of the Euler system in Example 15.11 The following two examples 
illustrate our observation that Richness 7^- rank(iV) = 0. In Example 15. 91 we have a rich system with 
rank(iV) = 1 with non trivial general solution depending on 1 arbitrary constant and 1 arbitrary 
function of 1 variable. In Example 15. 101 we have a rich system with rank(iV) = 2 and therefore only 
trivial solutions. 

Example 5.9. n = 3, rich, rank(iV) = 1, non-trivial. Consider a system with eigenvectors: 

i?i = (l,0,u 2 ) T , R 2 = (0,l,u 1 ) T , Rs = (0,0, -1) T . 

These vectors commute and therefore any system U.l\) with these as eigenvectors is rich. The 
only non-zero connection components are Tf 2 = r|i = — 1> such that the X-system implies a single 
algebraic relation: A 1 = A 2 . Thus this example illustrates the results in Section \4-3\ Taking into 
account this relation the differential equations become: 



nix 1 ) 


= diX^^dsX 1 


= 


ri(A 3 ) 


= <9iA 3 +-u 2 d 3 A 3 


= 


T2(X 1 ) 


= d 2 X 1 +u 1 d 3 X 1 


= 


r 2 (A 3 ) 


= 5 2 A 3 + u 1 9 3 A 3 


= 


^(A 1 ) 


-dsX 1 = 0. 





The general solution of this system is A 1 = A 2 = A, and A 3 = f(u s —u 1 ^), where X is an arbitrary 
constant and (p is an arbitrary function of one variable. 

Example 5.10. n = 3, rich, rankA^ = 2, trivial. Consider the vector fields 

Ri = [u\ u 2 , Of , R 2 = [-u 2 , u\ Of , R 3 = [-u 2 , u\ if . 

We are then searching for systems lil.l]) whose eigencurves are horizontal lines through the u s -axis, 
horizontal circles, and vertical helices. A calculation shows that these determine a rich frame. 
We can rewrite this problem in terms of Riemann invariants using the change of coordinates: 
p : M 3 — ► M 3 given by 

(5.23) p^ir^,^ = (u l ,u 2 ,u 3 ) = (rcos(9 + h),rsm(0 + h),h) . 

A calculation show that = —r and Z 2 3 = - . R follows from and Proposition \2.1\ that the 

only solutions are the trivial solutions A 1 = A 2 = A 3 = const. 
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